Vector-type four-quark interaction and its impact on QCD phase structure 
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Effects of the vector-type four-quark interaction on QCD phase structure are investigated in the imaginary 
chemical potential (/i) region, by using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with 
the extended Za symmetry. In the course to this end, we clarify analytically the Roberge- Weiss periodicity 
and symmetry properties of various quantities under the existence of a vector-type four-quark interaction. In 
the imaginary p region, the chiral condensate and the quark number density are sensitive to the strength of the 
interaction. Based on this result, we propose a possibility to determine the strength of the vector-type interaction, 
which largely affects QCD phase structure in the real /i region, by comparing the results of lattice simulations 
and effective model calculations in the imaginary fj, region. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

The progress in computer power has made it feasible to do 
realistic simulations in the lattice QCD for finite temperature 
(T) system without quark chemical potential (/i) [1]. As for 
p 2 > 0, however, lattice QCD has the well-known sign prob- 
lem, and then the results are still far from perfection; for ex- 
ample, see Ref. |2] and references therein. Several approaches 
have been proposed to solve the sign problem. One of them is 
the use of imaginary chemical potential, since the fermionic 
determinant that appears in the Euclidean partition function is 
real in the case; for example, see Refs. fUEll!!] and references 
therein. 

When physical quantities are available with lattice QCD in 
the imaginary // region, in principle it is possible to extrap- 
olate them to real fx, until there appears a discontinuity. Ac- 
tually, such an extrapolation was made for the phase transi- 
tion curve by assuming some analytic functions for the curve 
JUS]. This direct extrapolation may work for small real p/T, 
but its accuracy is quite unknown for large real fi/T [g]. This 
problem may be circumvented by the effective theory that can 
evaluate the partition function in both the real and imaginary // 
regions and reproduce the results of lattice QCD in the imagi- 
nary /i region, if such an effective theory is found. 

In the region of imaginary chemical potential fj, = iT9, 
Roberge and Weiss (RW) found [7] that the partition function 
Z(9) of SU(7V) gauge theory is a periodic function of 9 with 
a period 2ir/N, showing that Z(9 + 2nk/N) is reduced to 
Z{&) with the Zjv transformation for any integer k. The RW 
periodicity means that Z{&) is invariant under the combination 
of the Zjv transformation and another transformation 9^9 + 
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2nk/N, that is, under the extended Zat transformation |6] 

q -> Uq, A, -> UA v U~ l - l(d v U)U-\ 8^9 + ^,(1) 

where U (x, t) are elements of SU(A^) with the boundary con- 
dition U(x,f3) = exp(—2iirk/N)U(x,0), q is the fermion 
field, A v is the gauge field and (3 is the inverse of temper- 
ature T. Quantities invariant under the extended Z^ trans- 
formation, such as the thermodynamic potential (2(6) and the 
chiral condensate, keep the RW periodicity. Meanwhile, the 
Polyakov loop <P is not invariant under the transformation (fl}, 
since it is transformed as <P — > ^>exp(— i2irk/N), Such a 
non-invariant quantity does not have the periodicity, but this 
problem can be solved by the modified Polyakov loop defined 
later in Eq. (fT~2T > invariant under the extended Zat transfor- 
mation 1 6]. Roberge and Weiss also showed with perturba- 
tion that in the high T region d(2(9)/d8 and $(6) are dis- 
continuous at 9 — (2k + l)ir/N, and also found with the 
strongly coupled lattice theory that the discontinuities disap- 
pear in the low T region. The first-order RW phase tran- 
sition of d(2(6)/d6, that is, the quark number density n = 
—dQjdp, = i(d(2(9)/d9)/T was observed in lattice simula- 
tions nail. 

The Nambu-Jona-Lasinio (NJL) model H] is famous as 
a model for understanding chiral symmetry breaking, but it 
does not have the extended Z3 symmetry as well as the con- 
finement mechanism. In the previous work [6], we pointed 
out that the Polyakov-loop extended NJL (PNJL) model de- 
signed to have the confinement mechanism j p , UOl [ill Jl2l 

[3 53, LOl [H, LI3, LH, [M HH m, M M Ea \m 

possesses the extended Z3 symmetry. Hence the thermody- 
namic potential (2(9) of the PNJL model has the RW period- 
icity. Using the RW periodicity and the 6>-evenness of (2(6), 
(2(9) = (2(—9), we also found that #-odd (6>-even) quantities 
such as n (a) have first-order (second-order) phase transitions 
on lines 9 = (2k + 1)tt/3 in the 9-T plane. This result is 
consistent with the results of lattice QCD and even more in- 
formative. Thus, the extended Z3 symmetry is essential. The 
extended Z3 symmetry of finite 9 QCD is an extension of the 
Z3 symmetry of pure gauge QCD. 
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The confinement mechanism is a key to understand the 
physical real fj, world. Actually, it is predicted that the mech- 
anism largely shifts the critical endpoint 0271] to ward higher 
T and lower [i than the NJL model predicts IU7l HHH. In 
contrast, it is known 1I22I I24ll28ll29ll that the vector- type four- 
quark interaction (q"ffj,q) 2 largely moves the critical endpoint 
in the opposite direction, if it is newly added to the NJL and 
PNJL models. Thus, it is essential from the phenomenologi- 
cal point of view to determine the strength of the vector-type 
four-quark interaction. 



II. PNJL MODEL 



A. Formulation 



We consider the two-flavor PNJL model with the vector- 
type four-quark interaction. The PNJL Lagrangian is 



C =q(i^ v D v - m )q + G s [(qq) 2 + (qij 5 rq) 2 } 
-G v (q^q) 2 -U($[A],$n[A],T), 



(2) 



So far, the vector-type interaction was often ignored in the 
NJL and PNJL models. In the relativistic meson-nucleon the- 
ory lf30fl . meanwhile, the repulsive force mediated by vector 
mesons is essential to account for the saturation property of 
nuclear matter. Using the auxiliary field method, one can con- 
vert quark-quark interactions to meson-quark interactions; for 
example, see Refs. ll3ll L32L l33ll and references therein. In 
the hadron phase, furthermore, quarks have a large effective 
mass as a result of spontaneous chiral symmetry breaking, and 
then nucleons can be considered to be formed from such three 
heavy quarks, that is, three constituent quarks. It is then nat- 
ural to think that there exists a correspondence between the 
meson-nucleon interactions and the quark-quark interactions. 
In this point of view, it is quite likely that the vector-type four- 
quark interaction is not negligible and even significant in par- 
ticular at the finite quark-density region corresponding to the 
nuclear saturation density. 



In this paper we investigate roles of the vector-type four- 
quark interaction in the imaginary /i region and show that the 
chiral condensate u and the quark number density n are quan- 
tities sensitive to the strength of vector-type four-quark inter- 
action. We then propose to measure n as well as a with lattice 
QCD. If the strength of vector- type four-quark interaction is 
determined from the measured a and n, the PNJL model with 
the interaction can predict the phase diagram in the real \i re- 
gion with reasonable reliability. The previous proof [6] on the 
RW periodicity and the even/odd property of the extended Z3 
invariant quantities is not applicable to the case with vector- 
type four-quark interaction. We then prove the properties in a 
way different from the previous one ||6|]. 



In sec. m we present the PNJL model with the vector-type 
four-quark interaction, and derive the thermodynamic poten- 
tial invariant under the extended Z3 transformation. For the 
case of imaginary \l = iTO, it is proven that the extended 
Z3 invariant quantities exhibit the RW periodicity. For both 
of imaginary and real /i cases, their even/odd properties are 
determined. It is shown from the the RW periodicity and the 
even/odd properties that the RW phase transition appearing in 
9 = (2k + l)7r/3 is a family of first-order phase transitions in 
0-odd quantities and second-order ones in 9-even quantities. 
In sec. |III1 we show results of numerical calculations. Section 
HVl is devoted to Summary. 



where q denotes the quark field, mo does the current quark 
mass, t stands for the isospin matrix, and D v = d v — iA v — 
Jq/x for the complex chemical potential /1 = fj,R +iT9. Param- 
eters G s and G v represent the coupling constants of the scalar- 
and vector-type four-quark interactions, (qq) 2 and (qj^q) 2 , 
respectively. The Polyakov potential U, defined later in (O , 
is a function of the Polyakov loop $ and its Hermitian conju- 
gate <? H , 



1 

N' 



-TrL, 



with 



L(x) = Vexp 



<?H = ^TrLt, 



dr A4 (x, t) 



(3) 



(4) 



where A4 = iAo, N = 3, and V is the path ordering. In 
the PNJL model, <P and are treated as classical variables. 
The temporal component A4 is diagonal in the flavor space, 
because the color and the flavor space are completely sepa- 
rated out in the present case. In the Polyakov gauge, L can be 
written in a diagonal form in the color space lllll : 



diag(e l 



(5) 



where <f> a = <f> 3 + 8 /\/3, (j) b = -0 3 + fa/VH and <\> c = 
-(0a + <W = -20 8 /V3. 

The Polyakov loop <P is an exact order parameter of the 
spontaneous Z3 symmetry breaking in the pure gauge theory. 
Although the Z3 symmetry is not exact in the system with 
dynamical quarks, it still seems to be a good indicator of the 
deconfinement phase transition. Therefore, we use <P to define 
the deconfinement phase transition. 

Making the mean field approximation (MFA), one can get 
the Lagrangian density 



£mfa ^qii-f^D^ + j S v - (m + S a ))q 
-U(a,n)-U($,$n,T), 



(6) 



with 



a = (qq), n = (qj q), 

= "2 G v 71 . 



£« = ~2G f 



a. 



U = G s er - G v n\ (7) 
and obtain the thermodynamic potential fi by making the path 



3 



integral over the quark field 

d 3 p 



n = - 2N f 
+ im [i 



(2tt) 3 



3T(p) 



In [1 

+ U, 



3(<P + <Pa( 
3(^ H + $e 



-PE-(p)\ e -PE-(p) 



-3^-(p)i 



-3/3B+(p)i 



(8) 



with T(p) = v/p 2 + Af 2 , M = m + E s and £ ,± (p) = 
-E(p) ± jx. Here the effective chemical potential jl is defined 
by A = MR + - 2G v n. 
WeuseWofRef. 10]: 



U 



-f($n$) 2 , (9) 



' 2 T 7 



(10) 



Parameters of W are fitted to results of lattice simulations in 
the pure gauge system with finite T ll34l[35ll . In the pure gauge 
system the parameter To agrees with the critical temperature 



Td of deconfinement phase transition. Hence, as a reasonable 
choice, one can adjust To to Td = 270 MeV of pure gauge 
lattice simulations. However, the PNJL model with this value 
of To yields somewhat larger value of Td than ~ 173 MeV 
predicted by the full LQCD simulation S [H S. In this 
paper, we then take the rescaled value To = 190 MeV that 
gives T D =176 MeV in the PNJL calculation with fi = 0. 

As expected, fl is invariant under the extended Z3 transfor- 
mation, 



<P(0) -> $(0)e 



: 2nh 
1 3 



(11) 



Thus, the extended Z3 invariance of Q is held, even if the 
vector-type four-quark interaction is added to the PNJL Lan- 
grangian. In general, the symmetry persists, even if any sort 
of multi-quark interactions are added to the NJL sector. 

As mentioned above, it is essential to introduce the modi- 
fied Polyakov loop, 

V = e w $, *h = e~ l0 $H, (12) 
invariant under the the extended Z3 transformation ( fTTT >. 



J 



The extended Z3 transformation is then rewritten into 



,±i0 



V(9) $(0), $h(0) 



(13) 



and Q is also into 

n = -2N f 



d 3 p 



1 





3E{p) + ~ In [1 + We- f)E{p) e f) ^ 



0* t 1 



6 2 (T)T 4 



-/5B(p) e ~/3A„ 



-2/3£(p) 2/3/i„ 3iS 



+ e" 



-3 / 9S(p)„3/3Ar.„3ie 1 



3 ^ e -2^(p) e -2/3A„ e - 



■3i6 



e -3/3E(p) e -3/3A„ e -3ie 



3„3i9 



-(<V e 



-3i0\ 



-J7 



(14) 



where jl n = /ir — 2G v n, and the factor exp (±£30) is 
also invariant under the transformation ( TT3T >. Variables X = 
if', I^h, cr and n are determined from the stationary conditions 

dQ/dX = 0, (15) 

by solving the equations for X. The thermodynamic potential 
Q(9) at each 9 is obtained by inserting the solutions into dT4b . 
Such a calculation is possible for any complex \i. 

B. RW periodicity and even/odd property 

We begin with the case of imaginary chemical potential 
fi = iT0, that is, take /iR = in ( TBl i. The equations (TT~5T > 
depend on only through the factor exp (3i9), so that the so- 
lutions X are functions of the factor: 



Since 0(9) is obtained by inserting the solutions X into ([8]), 
£2(6) is also a function of exp (3i9), Hence, we get 

f2(fi+^) = f2(6), X(9+^)=X(9). (17) 

Thus, quantities invariant with respect to the extended Z3 
transformation have the RW periodicity. It is found from the 
proof shown above that the property is not influenced by the 
presence/absence of the four-quark vector-type interaction. 
This will be confirmed later through Fig[T] 

Next we discuss their more detailed symmetry properties. 
Taking the complex conjugate to (fT4l . one can find that the 
complex conjugate Q* has the same form as the original one 
]?, if a*, n*, and tjg are replaced by a and — n, <?h 
and if - , respectively. This indicates that the solutions X* of 
dQ*/dX* = are related to those X of dfl/dX = as 



X=X(e 3ie ). 



(16) 



a* = a, n* = -n, &* = tf H , = 



(18) 
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The last two equations represent that !fn is the complex con- 
jugate to !f , so we simply use W* instead of <P R m the case of 
imaginary [i. Meanwhile, the first two equations indicate 

d(2* _ dn dn* _ dn 

dnio drriQ ' dd dO 

because a = dQ/drriQ and n = —dQ/dpt = ij3d(2/dd. The 
relations ( fT9l ) are satisfied, only when (2 is real. Thus, (2 is 
real, so that n is pure imaginary and a is real. 

As shown in ( [To*l i, n is a function of cos(30) and i sin(30). 
This shows, together with the fact that n is pure imaginary, 
that n is 0-odd (odd under the transformation — > —6), since 

-n(cos(30),isin(30)) = n(cos(30), % sin(30))* 

= rt(cos(30), — i sin(30)) = rt(cos(— 30), i sin(— 30)). 

(20) 

Similarly, a is real and a function of 008(30) and zsin(30). 
This leads to the fact that a is 0-even, because 

cr(cos(30),isin(30)) = cr(cos(30), i sin(30))* 

= cr(cos(30),-isin(30)) = cr(cos(-30), i sin(-30)). 

(21) 

The potential J?(0) depends on through #(0), !f(0)*, 
cr(0), n(0) and e 3lS . We then denote 12(0) by (2(6) = 
(2(W(9),W(9)* ,n(9),e 3z9 ), where cr(0) is suppressed since 
it is 0-even and does not make any influence on discussions 
shown below. Equation ( fT4b keeps the same form under the 
transformation — > —0, if <f (0) and <f (0)* are replaced by 
>f (— 0)* and <f (— 0), respectively. This indicates that 

<f (-0) = if (0)* and if (-0)* = <f (0). (22) 

Using these properties and the fact that (2 is real, one can show 
that 

f}{6) =(f2(0))* = (2(&(6)*,&(6),-n(6),e- 3ie ) 

=(2(&(-6),&(-6)*,n(-6),e- m ) = (2(-6). (23) 

Therefore, we can summarize that (2 is a periodic even func- 
tion of with a period 2tt/3, (2(6) = (2(6 + 27rfc/3) = 
(2(—6). The chiral condensate a(6) is also a periodic even 
function of 0, while the quark number density n is a periodic 
odd function of 0. These properties are not changed by the 
presence/absence of the vector-type four-quark interaction, as 
understood from the proof shown above. 

The real (imaginary) part of <f is 0-even (0-odd), because 

Re[&(6)} = (lf(0)+#(0)*)/2 = Re[<f(-0)], 
lm[«f(0)] = (&(6)-&(6)*)/(2i) = -Im[W(-6)],(24) 

where use has been made of d22l . Thus, the real (imag- 
inary) part of !f is a periodic even (odd) function of 0. 
Similarly, the absolute value |!f | (phase (p) of the Polyakov 
loop is a periodic even (odd) function of 0, because |<f | = 
y/{Re[&]) 2 + (Im[!f]) 2 (0 = arctan(Im[<f]/Re[<f])). Obvi- 
ously, these properties persist also in the case of G v = 0. 



Since all quantities of our interest have the RW periodic- 
ity, here we consider a period < < 27r/3 in the 9-T 
plane. In the region, periodic even functions such as (2(9), 
cr(0), Re[!f (0)] and |!f | are symmetric with respect to a line 
= 7r/3. This indicates that such an even function X e has a 
cusp at = 7r/3, if the gradient dX c /d9\ g=7r / 3 ± e is not zero, 
where e is a positive infinitesimal. The phase transition ap- 
pearing in X c at = 7r/3 is of the second order. Meanwhile, 
Im[<f (0)], cf) and n are periodic odd functions of 0. This leads 
to the fact that these are discontinuous at = 7r/3, if the 
odd functions are not zero there. Thus, the phase transitions 
appearing in the odd functions at = tt/3 are of the first 
order. These are seen in the high T region, as shown later 
with numerical calculations. The phase transition appearing 
at = 7r/3 is called the RW transition. Thus, the RW phase 
transition is a family of first order transitions in 0-odd quan- 
tities and second order ones in 0-even quantities. This result 
is not changed by the presence/absence of vector-type four 
quark interaction. 

As mentioned above, (2 of ( TT4T > is invariant under the ex- 
tended Z3 transformation for any complex chemical potential 
M = MR + iT6. This implies that X has the RW periodic- 
ity, X(fi R + iT9) = X(fi R + iT(6 + 27rfc/3)). Actually, (2 
depends on only through the factor exp(3i0), even if /zr is 
nonzero. Hence, the solutions X of dfijdX — are func- 
tions of cxp(3i0). Thus, the concept of the extended Z3 trans- 
formation is useful for any complex chemical potential. 

Finally, we briefly consider the case of real n by taking 
= in (O; here note that «f = ^ and <f H = ^H- The 
complex conjugate potential (2* has the same form as (2, if a*, 
n* , If* and if^ are replaced by a, n, if and ifn, respectively. 
This indicates that the solutions X* of d(2* jdX* — are 
related to the solutions X of 8(2/ dX = as 

a* = a, n* = n, «f * = <f , >f h = ^h- (25) 

Thus, the four quantities are real. Furthermore, the first and 
second equations of d25l l are reduced to 

d(2* _ d(2 d(2* _ d(2 
dmo dmo ' dfj, R dfi R ' 

respectively. This indicates that (2 is also real. 

Equation ( TT4T > keeps the same form under the transforma- 
tion — > — /ir, if >f, 'fjj' cr and n are replaced by 
if, er and — n, respectively. This indicates that the solutions 
X(— /ir) of d(2(~HYi)/dX = are related to the solutions 
X(fx R ) ofd(2(p R )/dX = Oas 

•M-mr) = ^(Mr), ^(-Mr) = ^h(Mr), 
c(-Mr) = c(mr)) "(-Mr) = -n(^n). (27) 

Using < f27b in ( TT4l i. we can easily confirm that (2(fi R ) ~ 
(2(—fi R ), Thus, (2 and a are even functions of /1 in both the 
real and imaginary /i regions, so the chiral transition curve 
shown later is plotted in the /i 2 -T plane. 
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III. NUMERICAL RESULTS 

We proceed to numerical calculations. Since the PNJL 
model are nonrenormalizable, it is then needed to introduce 
a cutoff in the momentum integration. Here we take a three- 
dimensional momentum cutoff A. Hence, the present model 
has four parameters too, A, G s , G v in the NJL sector. Follow- 
ing Ref. JH, we take m = 5.5 MeV, A = 0.6315 GeV and 
G s = 5.498 GeV 2 that reproduce the pion decay constant 
/„ = 93.3 MeV and the pion mass M n = 138 MeV. That 
is, we pin down parameters other than G v at [i = 0, then we 
treat G v as a free parameter and vary its value in the range 
< G v < G s . 

Figure [TJ represents the RW periodicity of various quanti- 
ties calculated with/without the vector-type interaction. These 
graphs clearly shows that the RW periodicity is not affected 
by the inclusion of the vector-type interaction as argued ana- 
lytically. At T = 170 MeV, which is below the endpoint of 
the RW phase transition (point A of Fig.O, all quantities are 
smooth as functions of 8, whereas at T = 250 MeV, which 
is above the endpoint, the 0-even quantities, a and Re [<?"], ex- 
hibit second-order and the #-odd quantities, Im[n] and Im[<?], 
do first-order phase transitions. We note here that the present 
endpoint A is lower than in our previous works 1 6] because Tq 
in the Polyakov potential was rescaled in the present calcula- 
tion. 

Next we pay attention to the impact of the vector-type in- 
teraction on various quantities. Figure [2] represents the four 
quantities same as in Fig. [TJ for half a RW period, < 9 < 
7r/3, calculated with various strengths G v . Note that here we 
present the results at T = 300 MeV since the effect of the 
vector-type interaction is more conspicuous at higher temper- 
atures. First of all, the chiral condensate shown in Fig. [2f a), 
in particular around the RW transition at 6 = tt/3, is sen- 
sitive to G v . The calculated result means that the vector-type 
interaction reduces the constituent quark mass M. The vector- 
type interaction suppresses the effective chemical potential 
/I = iOT — 2iG v Im[n]. This directly affects Im[n], namely, 
Im[n] decreases as G v increases. Although the indirect effect 
through M, mentioned just above, works oppositely, the direct 
effect is stronger and therefore survives. Consequently effects 
of the vector-type interaction on Im[n] shown in Fig.|2|b) as 
well as on a are significant. Figures |2J C ) an d (d) indicate that 
the effects on the modified Polyakov loop are negligible be- 
cause it is fully determined by the Polyakov potential. The 
physical value of G v is still an open matter although its effect 
is significant in the physical real-/i phase diagram (see Fig. [4] 
below). Since the present calculation shows that its impacts 
on the chiral condensate and the quark number density are ap- 
preciable in the imaginary-/! region where both lattice simula- 
tions and effective model calculations are available, we expect 
that it would be possible to determine a physical strength of 
G v by comparing both the results, although the former is not 
available yet for the case of the two flavor QCD. 

We report the impact of the vector-type four-quark inter- 
action on the n 2 -T phase diagram ranging from negative to 
positive /i 2 in Fig. [3] since the order parameters a are func- 
tions of /i 2 as proved analytically above. This figure indi- 
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(6-odd quantities) are for T = 250 MeV, while smooth ones are for 
T = 170 MeV. 
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Fig. 2: The impact of the vector-type four-quark interaction on the 
four quantities same as in Fig.[j]at T = 300 Me V for a half Roberge- 
Weiss period. 



cates that effect of the vector-type interaction on the phase 
diagram is visible in the real-/i part; it has been known that 
the critical endpoint (CEP; point C) moves to lower T and 
higher /j, as G v increases and eventually disappears at G v « 
0.38G S JHHimiH]. In contrast, its effects are less in the 
imaginary-/i part of the phase diagram although the values of 
a and lm[n] are sensitive to it. This can be understood from 
the fact that the effect of finite G v on fi is important when T 
is low whereas the RW transition occurs at higher T. We note 
here that this figure graphs the chiral transition curve although 
the deconfinement transition temperature in the imaginary-/! 
is lower. We confirmed that their difference decreases as 9 
increases due to finite G v . 
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Fig. 3: The phase diagram in the /i 2 -T plane, (a) for G v = and 
(b) for G v = 0.25G S . The point A is the endpoint of the RW phase 
transition, B is the crossing point of the ordinary chiral crossover 
transition curve and the RW phase transition line, C is the chiral crit- 
ical endpoint, and D is the chiral transition point at T — 0. 

Finally, in order to look into the correlation between the key 
quantities in the imaginary- and real-/i regions, we present 
in Fig. [4] the relation between the chiral condensate at 9 = 
7r/3 — e and T = 250 MeV in the imaginary-/! region, which 
can be directly measured in lattice simulations, and the posi- 
tion of the chiral CEP in the real-/! region, which cannot be 
measured. Both of them are functions of G v ; here e is a pos- 
itive infinitesimal. The larger G v is adopted, the smaller —a 
becomes. The right-most point corresponds to G v = while 
the left-most does to G v — 0.30G S . Since the correlation 
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shown in Fig. [4] is evident, we hope this serves to determine 
the physical value of G v . 
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Fig. 4: The correlation between the chiral condensate at 9 — 7r/3 — e 
in the imaginary-^i region and the position of the chiral critical end- 
point (CEP) in the real-/i region, where e is a positive infinitesimal. 
These curves correspond to < G v < 0.30G S . 



Nambu-Jona-Lasinio (PNJL) model possesses the symmetry, 
even if the vector-type four-quark interaction is newly added 
to the NJL sector of the PNJL Lagrangian. Thus, PNJL is a 
suitable model to analyze not only the real /i region but also 
the imaginary one. 

Using PNJL, we have analyzed the effect of the vector-type 
four-quark interaction in the imaginary /i region and investi- 
gated the correlation between the position of the critical end- 
point (CEP) in the real /i region and several quantities in the 
imaginary /i region. In the real /i region, the position of CEP 
is known to be sensitive to the strength of the vector-type four- 
quark interaction. Meanwhile, the present analysis shows that 
in the imaginary /i region the vector-type four-quark interac- 
tion largely changes values of the chiral condensate and the 
quark-number density. This indicates that these quantities in 
the imaginary /i region are suitable to determine the strength 
and therefore the QCD phase diagram in the real /i region. 
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